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We point out that in inclusive B — > Xs 



decay an angular decomposition provides a third 



{q dependent) observable sensitive to a different combination of Wilson coefficients than the rate 
and the forward-backward asymmetry. Since a precise measurement of dependence requires large 
data sets, it is important to consider the data integrated over regions of . We develop a strategy 
to extract all measurable Wilson coefBcients in B ^ Xst^l~ from a few simple integrated rates 
in the low region. A similar decomposition in B ^ K*l^t~ , together with the B — > K*'y rate, 
also provides a determination of the Wilson coefBcients, without reliance on form factor models and 
without having to measure the zero of the forward-backward asymmetry. 



I. INTRODUCTION 

At the scale of B meson decays, flavor changing inter- 
actions at the electroweak scale and above are encoded 
in Wilson coefflcients of operators of dimension five and 
higher. The main goal of the B physics program is to 
make overconstraining measurements of the magnitudes 
and phases of these coefficients and thereby search 
for deviations from the standard model (SM). 

The b st^£~ process has been observed both 
in inclusive B Xs£^i~ ^ |^ and exclusive B 
K(*)e+l- g I decays. Throughout this paper we as- 
sume the SM except where explicitly stated otherwise. 
We also neglect the strange-quark and lepton masses. 
Two observables that have been extensively discussed for 
inclusive B Xsi^i~ decay are the spectrum j^] and 
the forward-backward asymmetry (or, equivalcntly, 
the energy asymmetry H). At lowest order, they are 
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Here q^ = [pp 
s = q^/ml, z = 



+ is the dilepton invariant mass, 

cobO, and 
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In B^ or B^ [B^ or _B+] decay, 9 is the angle between 
the ^+ [t^] and the B meson three- momenta in the t^t'^ 
center-of-mass frame. The Wilson coefficients Cy^g.io 
contain short distance information. Beyond tree level 
they are effectively dependent and complex, and re- 
ceive different contributions in inclusive and exclusive de- 
cays, as will be discussed below. If there were very precise 



data on _B ^ Xst'^i~ , one could extract the individual 
Wilson coefficients from the q^ dependence of dV /dq^ and 
the zero of the forward-backward asymmetry. As long 
as the measurements are limited by experimental uncer- 
tainties, it is important to find the most effective ways to 
extract the short distance information from a few simple 
observables integrated over q^, accessible with a limited 
amount of data. 

In Section || we discuss general aspects of an an- 
gular decomposition, which gives three observables, 
Ht,a,l{<I^)- To extract short distance information from 
these, we separate out the part of the Wilson coefficients 
sensitive to new physics in a q^ and /i independent man- 
ner, and propose to compare measurements of these with 
their SM predictions. Section [II investigates inclusive 
B Xsi'^i^ decay, while Section [rv| d eals with the ex- 
clusive B K*i'^i~ mode. Section |V| contains our con- 
clusions. Many analytical results and numerical inputs 
are collected in the Appendices. 



II. ANGULAR DECOMPOSITION AND 
DEPENDENCE ON WILSON COEFFICIENTS 

The double differential decay rate in q^ and z = cos t 
for either inclusive B Xs£~^i~ or exclusive B — 
K*i~^i~ decays can be written as 



dg2 dz 



-[{l + Z^)HT{q^) + 2zHA{q^ 
+ 2{1 -z')HL{q')]. 



(4) 



The functions Hi{q^) are defined to contain the full q^ 
dependence of the rate and are independent of z. They 
can be extracted from the d?T /dq^dz distribution either 
by a direct fit to the z dependence or by taking integrals 
of z. As special cases we have 
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For Hl the hadronic current is longitudinally polar- 



ized, so the rate goes like sin 9 — 1 



For Hx and 



Ha the hadronic current is transversely polarized, with 
Ht containing the contributions from purely vector and 
axial-vector Icptonic currents and Ha containing the in- 
terference between vector and axial-vector leptonic cur- 
rents. In the combinations Ht ± Ha the hadronic and 
leptonic currents have the same (opposite) helicity, giv- 
ing the usual (1 ± cos6')'^ = (1 ± z)^ dependence. This 
decomposition is a common tool in the analysis of exclu- 
sive semileptonic decays (e.g., B D*tv, piv, K*i~^£^) 
and of decays to two vector mesons (e.g., B <pK*, 
J/ijjK*). In analyzing inclusive B Xgi^i^, it should 
not be harder than measuring Afb- 

We introduce a scheme to separate certain SM contri- 
butions to the rate from terms that are most sensitive to 
new physics. We define modified Wilson coefficients 



Ct = Crip) [rnb{^I)/ml^] 

Cg = Cg(/l) + . . . , 



(6) 



The ellipses denote a minimal set of perturbative correc- 
tions, such that C7.9 are ^ independent and real in the 
SM. These coefficients are given explicitly in Eq. (A2) 
in Appendix ^ to 0{as). The decay rate also depends 
on SM contributions that are not contained in Cr^g.io- 
These are discussed in Sees. Ill and [V, Of these, the 



dominant contributions are from the four-quark opera- 
tors Oi,2, which are expected to be given by the SM. We 
regard the C7, Cg, and Cio as the unknown parameters 
that need to be extracted from experimental data, and 
compared with the SM or new physics predictions. 

We do not expand the MS &-quark mass, Wibifi), in 
C7, because it always comes together with O7 and they 
should be renormalized together. This also makes the 
perturbative expansions better behaved. In addition, as 
indicated in Eq. (||), we use the IS* scheme for all other 
factors of rrib (and rric as well) , which also improves the 
perturbative expansions. We will drop the superscript 
IS hereafter when the distinction is unimportant, but 
use ml^ for mt everywhere (except nib, of course). 

At subleading orders in as and l/wf,, the dependence 
of the rates on the Ci's will be different in inclusive and 
exclusive decays. To simplify the explanation of our main 
points, in the remainder of this section we neglect the 
contributions from operators other than 07,9.10. Then, 
at leading order, the i?i's defined in Eq. (|^) have the 
general structure 
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HA{q^) cx -4(1 - s fsCioiCg + ^ C7) 
HL{q^)^{l-sf[{Cg + 2Crf+Cfo] 



(7) 



Comparing Eqs. (|i|) and (|^) shows that splitting dT/dq^ 
into HL{q^) and Hxiq^) separates the contributions with 



different q^ dependences, providing a third independent 
observable, which has not been studied so far in inclu- 
sive B -> Xsl+l~. For exclusive B -> K*£+e~- it is 
equivalent to the fraction of longitudinal polarization 
Fl = Hl/{Ht + Hl) measured by Babar 

If one does not resolve the details of the hadronic sys- 
tem, neglects effects proportional to mj/ml, and does 
not measure the lepton polarization (which might be ac- 
cessible in S — > XsT'^T" |T^, though this is challenging 
at best), then the only observable linear combinations of 
the Wilson coefficients are those given in Eq. (R). This 
is necessarily the case for inclusive B Xse~e^ and 
B ^ Xsp,~^ fi~ decays. 

In exclusive B K*£+£~ decay, if the K* Ktt 
decay is reconstructed, there are two additional observ- 
able angles. These are 9k* , which is the analog of 9 for 
the Kt: system, and 0, which is the angle between the 
Kit and the planes (using the notation of Babar [|j , 
which follows Ref. [0). Once one integrates over 0, even 
if the 9 and 9k-' distributions are not integrated over, 
the rate depends only on the three linear combinations 
in Eq. (^. Keeping the dependence would give rise to 
two further linear combinations , and it would require 
a more detailed study to test whether the measurement 
could benefit from not integrating over cj). 

As mentioned above, instead of relying on the full q^ 
dependence, we want to integrate over as large regions of 
q^ as possible to extract the Wilson coefficients from the 
simple integrals 



H^iql ql) 
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We restrict our discussion to the low q region, 1 GeV < 
ql,q2 < 6GeV^, since it is theoretically clean and con- 
tains a large part of the rate. The interference of the 
J/^ contribution with the short distance rate is a sig- 
nificant contamination at higher values of q^, while the 
rate for q^ > m^, « 14.2 GeV^ is significantly smaller. 
Ultimately, the measured tail of the long distance contri- 
bution will determine the optimal upper cut on q^. 

For HL{q^) the hadronic current is longitudinally po- 
larized, so the C7 contribution is not enhanced by a 1/s 
pole, and is numerically small. Since the Wilson coeffi- 
cients in Hl combine into a q^ independent overall fac- 
tor, there is no gain in considering the (f' dependence of 
HL{q^). Thus, to get maximal statistics one should use 



i/L(l,6) cxC2 + 
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4C7Cg -f 4C^, 
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which is dominated by the C| + C^q term. 

In contrast, in Hxiq^) the different contributions have 
a hierarchical q^ dependence, i.e. 



Hriq^) cx -C 
s 



+ AC7C9 + s{Cl+Cf^) 



(10) 



In this case, integrating H-riq^) over q^, the C7Cg interfer- 
ence term in Ht{1, 6) is as important as the Cq+Ciq con- 
tribution, because the latter vanishes for q^ — s- 0. Thus, 
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with enough statistics, it will be worth splitting Ht{1, 6) 
into two integrals, Hxil^qf) and HT{qf,6), giving ac- 
cess to two independent combinations of Wilson coeffi- 
cients. The value of qf should be chosen to give compa- 
rable statistics for the two integrated rates. We will use 
Ht(1,3.5) and iJT(3.5,6) below. 

Finally, _ff^(l,6) oc ApB has comparable C7C10 and 
C9C10 contributions. To separate them, one can split 
^^(1,6) into HA{l,qi) and HA{qj,Q) similarly to Ht 
above, where the precise value of q^ can again vary. We 
will use 3.5) and i?^(3.5, 6), which give reasonably 

independent linear combinations of C7C10 and C9C10 . We 
do not normalize iJ^(q^) by dF/dg^, as is often done for 
dA-p-e/dq^. This would suppress the low q^ region, with- 
out providing any real advantage. If necessary, HA{q^) 
could be normalized to the rate in a certain q^ window. 

Usually the zero of HA{q^) oc dApB/dq^, which occurs 
in the vicinity of — —2mf^{C-j / Cq) , is advocated to 
determine Cj/Cq. As discussed in the introduction, mea- 
suring requires very large data sets. Measuring two 
integrals of Ha (q^ ) as described above may be a simpler 
way to achieve similar sensitivity with less data. 



III. INCLUSIVE B 



In this section we consider the inclusive decay B — > 
Xgi'^i^ , working to what is usually referred to as next- 
to-next-to-leading order (NNLO). We define the effective 
coefficients 



cr'('z') = c9 



F,{q^) 
F,{q^) 



-G7{q') 
-G,{q^) 



(11) 



such that all terms on the right-hand side of Eq. ( pi] ) are 
separately /i independent to the order we are working at. 

The functions Fi{q^) and Gi{q^) are calculated in 
the SM. The Fr^g{q'^) contain perturbative contributions 
from the four-quark operators Oi_6 and the chromomag- 
netic penguin operator, Og, while the GT^g{q^) contain 
nonperturbative Oil/m^) corrections involving the four- 
quark operators |12| . The latter can be included in a sim- 
ple form for any differential rate, but the final results have 
to be reexpanded so that 0{as/m'^, l/m*) terms are not 
kept. The explicit expressions are given in Appendix 
In the small q^ region (well below the cc threshold), the 
Cy"g'(g^) have only small imaginary parts and modest q^ 
dependences, which arise only from Oi-6,8 and are fully 
contained in Fi{q^) and Gi{q^). Therefore, all Ci are real 
numbers in the SM, which has the advantage of reducing 
the number of parameters to three. 

At NNLO, we include the corrections to the Fi{q^ 
up to 0{as) where they are known analytically 
p^, |l^, and to the Wilson coefficients of Oi_6,8 enter- 

191 [20[ |. At leading order in 
e Fg{g^) is nonvanishing. We 
also include the 0{as) corrections to the matrix ele- 
ments of O7 9jn |2^ and the nonperturbative 
0{l/ml) - 



M corrections. The short distance contri- 
butions to the Hi{q^) can be written as [s = q^ /{ml^Y] 



ing the F.(g2) 0, |T8| 
Q!s, F'j{q^) vanishes, whr 



Hriq^) = 2Fo(m^^)^ (1 - sfs^C^^']^ + Cl^)hf{s) + l|Crf + - Re(Cr'*Cr')4'(s) 
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-4Fo(mi^)3(l-s)2sCioRe 
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HL{q') - Fo(mi^)3 (1 - s)' {\Ci^'^'\' + C!„)hf{s) + 4\Cr'\'hY {s) + iRc{a 
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The functions /i^(s) are defined to have only a residual s dependence entering at higher orders in and l/rrii,, i.e. 

1 



2^-^(.) 



(13) 



with Cf — 4/3. The ujI{s) and xli^) containing the 0{as) and 0{\/m1) corrections are given in Appendix]^. We 
neglect the finite bremsstrahlung corrections due to four-quark operators, i?j^'°™''(g^), which are not known for the 
double differential rate d^F/dg^dz. This should be a safe approximation, since in the known cases of dF/dg^ |2^ and 
d^FB / dg^ ||2^ they are at or below the one percent level. 

To make numerical predictions we use the input values collected in Appendix 
to better than one percent, there would not be a significant benefit relative to other uncertainties in normalizing the 
rates to F(B ^ Xiv). We find 



Since the value of is known 



i/T(l,3.5)/Fo = 37.18 (C^+Ci^ 



15550. + 1409. C7C9 - II9.4C9 - 2795. C7 ■ 



69.65. 



i/T (3.5, 6)/Fo = 59.76 {Cl + C^o) + 1067. C7C9 + 5008. - 76.25 Cg - 778.2 C7 - 67.72 , 
Ha{1, 3.5)/Fo = -Cio (7O.I9C9 + 1401. C7 - 121.5) , 
ffA(3.5, 6)/Fo = -Cio (III.8C9 + 1051. C7 - 80.91) , 
Hl{1, 6)/Fo = 315.2 {Cl + Cjo) + 1377. + 1299. C7C9 + 33.41 Cg -I- 86.86 C7 - 72.87 . 



(14) 
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FIG. 1: (color online) Constraints in the Cg — Cio plane. Left: r(l, 3.5) and r(3.5, 6) [light and medium (gray) annuli], HA{i, 3.5) 
and Ha{3.5, 6) [light and dark (green) regions bounded by hyperbolae]. Right: Ht{1, 6) [dark (blue) annulus], Hl{1, 6) [medium 
(orange) annulus], //a(1,6) [light (green) region bounded by hyperbolae], and for comparison r(l,6) [light (gray) annulus]. 
The black dots show the assumed (SM) central values. The combined constraints are also shown [small light (yellow) regions]. 



The major uncertainties in Eqs. (|TJ) arise from higher 
order perturbative corrections, TOf, and rric- Varying the 
renormaUzation scale between mb/2 and 2mb, we get less 
than 5% uncertainty in the coefficients of the dominant 
terms in Eq. . Since the difference mt — rric is known 
precisely [^, we vary rrib and rric in a correlated manner, 
which gives a 1-5% uncertainty. The uncertainties from 
other input parameters and higher order corrections in 
l/mj, are much smaller. The uncertainties from the elec- 
troweak matching scale, /ip, and the top-quark mass, mj, 
in Eq. ( p^ ) are negligible, because they primarily enter 
via the values of the d. Using Eq. (|lj) and the values of 
Ci from Appendix we obtain the SM prediction for the 
B X,i+i- branching ratio for 1 GeV^ < < 6 GeV^ 



rBr(l, 6) = (1.575 ± 0.067[^] ± 0.051[™,,„^] 
± 0.041 [,„,]±0.019[„,j) X 10-^ 



(15) 



This agrees well with Refs. |20| and with Ref. [|8|, 
which also uses the IS* scheme. 

To illustrate the improvement in determining the Wil- 
son coefficients one can obtain by separating T into Ht 
and Hl, we scale the current measurements |^ to 
1 ab""'^ luminosity, which gives about 10% statistical un- 
certainty for r(l,6). We assume that the Hi are mea- 
sured with the central values given by the SM. The sta- 
tistical error of Ht and Hl is obtained by scaling by 
the number of events compared to r(l, 6). In the case of 
Ha we take the same absolute statistical error for 3 / AH a 
as for the total rate integrated over the same g^-region. 
The reason is that ?>/AHa corresponds to the difference 
between the rates for positive and negative cos 9, which 



has the same absolute statistical error as the sum. To 
this we add in quadrature a 20% systematic uncertainty 
for all Hi, to account for experimental systematics and 
theoretical uncertainties. From each observable's total 
error we build for the individual and combined con- 
straints, and Figs. |l| and ^ show the Ax^ = 1 regions 
in the Cg — Cio plane. Since B will always be 

measured with higher precision than B — > Xst^t-~ , we 
consider the value of Cj to be known from B — > Xs^ and 
assume its sign is negative as in the SM (since there is 
an overall sign ambiguity). 

On the left-hand side in Fig. |l| wc show the constraints 
from r(l, 3.5), r(3.5, 6) [hght and medium (gray) annuli], 
Ha{^, 3.5), and iJyi(3.5, 6) [light and dark (green) regions 
bounded by hyperbolae]. This plot shows that splitting 
r(l, 6) into two regions is not really useful because a very 
similar linear combination of Wilson coefficients is con- 
strained. As is well known, splitting 7f^(l,6) into two 
regions is very useful, since different combinations of co- 
efficients are constrained by each region. (This is also 
the reason why the zero of the forward-backward asym- 
metry is interesting to study.) The plot on the right in 
Fig. |l| shows that splitting r(l,6) into Ht{^,&) [dark 
(blue) annulus] and Hl{\,&) [medium (orange) annulus] 
gives a very powerful constraint. This shows the power of 
separately measuring Ht and H^ as advocated in the in- 
troduction. The observables shown in this figure are suffi- 
cient to extract the absolute values \Ci \ and the sign of Cg 
relative to Cy. The constraint from r(l, 6) is also plotted 
[light (gray) annulus] , which shows that separating F into 
Ht and H^ significantly improves the constraints. How- 
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FIG. 2: (color online) Constraints in the Cg — Cio plane. Left: Ht{1,(3) [dark (blue) annulus], Hl{1,6) [medium (orange) 
annulus], Ha{1,3-5) and Ha{S-5,6) [light and dark (green) regions bounded by hyperbolae]. Right: as on the left, but Ht 
split into -ffr(l,3.5) and Ht{S.5,6) [dark and light (blue) annuli]. The right plot also shows the constraint from the two Ht 
observables alone [small black (dark blue) regions]. All other notations are as in Fig hi. 



ever, because of its large relative error, _ff^(l,6) [light 
(green) region bounded by hyperbolae] does not provide 
a good constraint. 

The left plot in Fig. |^ shows that spUtting Ha into 
two regions gives sensitivity to the sign of Ciq. Measur- 
ing Ha{1,3.5) and Ha{3.5,6) can distinguish between 
the positive and negative solutions for Ciq. Combined 
with the tighter constraints from spUtting F into Ht and 
Hl, spUtting Ha into these two regions gives informa- 
tion similar to the zero of ApB- The plot on the right 
in Fig. H shows that splitting 7Jt(1,6) into ffT(l,3.5) 
and iJT(3.5,6) can further overconstrain the determina- 
tion of the Wilson coefficients. The black (dark blue) 
region in the right plot in Fig. ^ shows the combined con- 
straint from only the two Ht integrals. The requirement 
that the Ht constraint overlaps with it effectively pro- 
vides a consistency test on the value of C7 extracted from 
B —> Xs7. Such overconstraining determinations of the 
Wilson coefficients provide model independent searches 
for physics beyond the SM, complementary to and possi- 
bly more effective than constraining specific new physics 
scenarios (e.g., Wilson coefficients with flipped signs, 
complex values, or adding right-handed operators). For 
example, if there are operators with right-handed helicity 
structure, they will affect Ht and Ht differently, because 
of the different polarizations. 

It would also be interesting to explore experimentally 
whether the influence of the J/ip resonance turns on at 
similar values in Ht, Ha, and Ht- Since we have 
no information on the J/^ polarization in inclusive B — > 
J/ ^pXs decay, it is possible that the upper cut on can 
be extended past 6 GeV^ in some (but maybe not all) 



of these observables, which may improve the statistical 
accuracy of the measurement. 



IV. EXCLUSIVE B K't+l' 

We now turn to the exclusive decay B K*i^£^. 
While the theoretical uncertainties are larger than in the 
inclusive analysis, measuring the exclusive mode is sim- 
pler and it may be the only possibility at LHCb. Com- 
pared with the inclusive decay, the exclusive measure- 



ments go closer to q 



8.1 GeV'' at Belle M and 



8.4 GeV^ at Babar §. In this region of phase space the 
energy of the K* varies only between 1.9 GeV < Ek' < 
2.7 GeV, which helps control some theoretical uncertain- 
ties. In our general discussion we will consider for sim- 
plicity 0.1 GeV^ < <7^ < 8 GeV^, where the precise value 
of neither limit is important (the lower limit can be re- 
placed by any experimentally appropriate value above 
Arnj). However, for comparisons of our results with the 
data, we use the limits used in the experimental analysis. 

In this section we explain that, similarly to the in- 
clusive decay, all the information obtainable can be ex- 
tracted from a few integrated rates. To obtain the most 
information on the ratios of Wilson coefficients. Belle jj] 
performed a maximum-likelihood fit to the double differ- 
ential distribution (PT/dq^dz. However, since the theo- 
retical predictions change for the double differential rate 
as they are being refined, we think that a few integrated 
rates will be very useful to compare the theory with the 
data. In addition, results from different experiments are 
more straightforward to combine for these partial rates. 
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In the heavy quark hmit, soft-colhnear effective theory 
(SCET) Hi relates the seven full QCD form factors that 
describe B K*£^£~ to fewer functions. We follow 
the notation of Ref. and denote these by || and 
II . The II obey the form factor relations js^, while 



d II violate them. Whether /( is ©(a^) or 0(1) is 
subject to discussion, and so at the present time the as 
corrections in exclusive decay are not fully known. 
The angular decomposition for B K*i'^£~ is 

MM 



cr' ais) + ^ ^ + (i - .) cm 

: niB 



HL{q^) ^ ^Fornix' [C 



cxcl '^h 



IS 



TUB 



(16) 



In this section s — /rn?g. In Eq. (|T^, p — mj^t/m^ ^ nonperturbative parameters, and neglect 0{as) correc- 
0.03 and A = [(1 - s)^ - 2p (1 + s) + p^]^/^. In analogy tions to C±\s), which are partially known. These same 



to Eq. (|l l|), we have defined 



cr'('?') = C7 + i^7(g') + 0(«.): 
CT'(9')=C9 + F9(g2) + 0(a,), 



(17) 



Note that there are additional 0{as) corrections beyond 
Ft, 9(1'^) [Hi 01 and we do not attempt to address power 
suppressed terms. 

Without lattice QCD (LQCD) input, model calcula- 
tions, or nonleptonic decay data to constrain Cy ^'i we 
cannot learn about the C^'s from Hi. However, the mag- 



nitudes of the form factors C 
the B K*j rate. 



G 



(J) 



can be constrained from 



TiB ^ K*j) ^ ^ ^ \V,,V;fm-U<'r (1 - P) 



Stt^ 47r 

x\cr\o)f [a(o) + cl( 



(18) 



In Eqs. wM, (16), and (03) we display the kinematical 



p dependences, but neglect the "dynamical" ones, which 
would, for example, multiply Cl(0) in Eq. (|l^ ) by {1 — p). 

The ( and ("^ form factors are functions of the K* 
energy, Ek' = (mB/2)(l — s -\- p). In the heavy quark 
limit, a convenient parameterization is 



Cf(0) 

(l-.)2 



[l + 0{as,A/EK')] , (19) 



i.e., the form factors' leading E dependence is 1/E'j^,. 
Since in the 0.1 GeV^ < q^ < 8.4 GeV^ region Ek' varies 
only over 1.9 GeV < Ek* < 2.7 GeV, we do not ex- 
pect large deviations from this limit. For example, C.± 
can have additional logarithmic dependence on Ek* , 
which is roughly constant over this small region. 

In the literature C'/C is often treated as 0{as). The 
zero-bin shows that this is not the case parametri- 
cally, so we treat both Cx(0) and C±{0) as independent 



considerations also imply that the zero of the forward- 
backward asymmetry in B ^ K*£^£^, does not 
necessarily provide as precise and as model independent 
a determination of C7 /Cg as claimed in much of the liter- 
ature. Moreover, even after 5 years of LHCb data taking 
(10fb~^), one expects cr^q^^) « 0.5 GeV^, which would 
determine C^/Cg only with an approximately 13% error 
in the SM |7). 

Some of the known 0{as) corrections to Eqs. ([l^ ) and 
( p^ ) can be included in Eq. (p7|). However, since we treat 
C±{0) and Ci(0) as independent unknowns and deter- 
mine them from the data, only the part of the as correc- 
tions that causes these form factors to deviate from their 
asymptotic s dependences in Eq. ( p^ ) will introduce er- 
rors. Since in the 0.1 GeV^ < q^ < 8GeV^ region that 
we concentrate on, < s < 0.3, we expect that the ne- 
glected 0{as) terms do not introduce a dominant error. 
Moreover, they could be added to our analysis. 

From the ratios of the four observables 



r(S^X*7), i/T(0.1,8), Ha (0.1, 4), i/A(4, 



one can extract 

Cg 
C7' 



c7 



cm 



a(o) + ci(o) 



(20) 



(21) 



If one does not want to use the B — > K*j data then 
ratios of Ht{0.1,4), Ht{4:,8), i7A(0.1,4), and i?A(4,8) 
could also be used. The ratio r is of great interest for 
heavy quark theory. One expects r to be roughly similar 
in size to the same ratio involving the Cy'^'' form factors, 
which enters the determination of the unitarity triangle 
angle 7 (or a) from charmless two-body B decays. 

At present there is insufficient published data to carry 
out this analysis to determine the Wilson coefficients. 
However, T{B — ^ K*^) is well measured and 
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FIG. 3: (color online) The ratio R{0.1, 8.4) defined in Eq. (H) 
as a function of r [dark (blue) curve] , assuming the SM values 
for all Ci. The medium (orange) curve includes the 0(1) cor- 
rection from Fglq^), and the light (green) curve in addition 
includes 0{as) corrections from ^7,9(5^) in Eq. (p^). The 
dashed [dotted] line shows the central value [la upper bound] 
from the Babar measurement [^]. 



polarization fraction in this region. 

Until precise unquenched LQCD calculations of the 
B K^*^ form factors for small become available, 
the method outlined above may provide the most accu- 
rate extraction of short distance information from B 
K*£~^£~. With more statistics in the future it should 
become possible to determine the quantities in Eq. (pl|), 
providing new tests of the SM and insights into the the- 
ory of hadronic B decays. 

We did not consider B Ki^£~ decay, because 
B Kj is forbidden by angular momentum conser- 
vation, so it is not possible to learn about the short 
distance physics from this mode without using a deter- 
mination of the corresponding form factors from lattice 
QCD or model calculations. (This is similar to the case 
of explained above.) To proceed by using 

model independent continuum methods, one would have 
to use B ttID data combined with SU (3) flavor sym- 
metry or the two-body charmless nonleptonic decay data 
to constrain the B K form factor. 



ffT(0. 1,8.4) can be obtained from the Babar measure- 
ment [^ . (Belle [^ does a maximum-likelihood fit to ex- 
tract information on the Wilson coefficients, so we cannot 
use their fit result projected on to the forward-backward 
asymmetry as a direct measurement of Ha in the corre- 
sponding bins.) Assuming the SM for Cy.g^io, we can use 
T{B -> K*j) and i/T(0.1,8.4) to constrain r. The ratio 



T{B -> K*-/) 



2 / 2 



ds 



A3 



C.2 



^(1 



C7 



(1 



(l_p)3(l_,)4 

r) + ^(1 - sr) 

S TTlB 



(22) 



as a function of r is shown in Fig. || for _R(0. 1,8.4). 
The dark (blue) curve shows Eq. ([2^), the medium (or- 
ange) one includes the leading corrections to C|™'(q^) 
from Fg{q'^), and the light (green) curve includes in addi- 
tion the 0{as) corrections from FT^g{q^). The significant 
change is dominantly due to the large ^7(0) contribution 
to Cif'"'\0), also observed for r{B K*-f) ||, M. This 
shows that a complete understanding of the as C^correc- 
tions to these exclusive decays is very important. Until 
this is achieved, a determination of Cq/Cj, C10/C7, and r 
without the r{B K*^) data, using only two bins of 
each Ht and 7?^, as mentioned after Eq. (^l|), may be 
theoretically cleaner. 

The central value, i?(0.1,8.4) = 1.55 x 10"^ ||, and 
the 1(7 upper bound are shown in Fig. || by dashed and 
dotted lines, respectively. Since the error is still large, at 
the moment one cannot make a statistically significant 
statement about the size of r. Nevertheless, we expect 
that (if the SM is valid) the central value of i?(0. 1,8.4) 
should go up, probably via an increase in the transverse 



V. CONCLUSIONS 

In this paper we pointed out that in inclusive B — > 
Xst'^£~ decay an angular decomposition provides a third 
q^ dependent observable in addition to the total rate and 
the forward-backward asymmetry. Splitting up the rate 
into transverse and longitudinal parts, proportional to 
1 -f cos^ 9 {Ht) and 1 — cos^ 9 (Hl), one gains access to 
a third independent linear combination of Wilson coeffi- 
cients. This requires measuring no additional kinematical 
variable besides q^ and cos6', which are already studied 
by Babar and Belle. Without doing more complicated 
analyses, it will improve the determination of the rele- 
vant Wilson coefficients and the sensitivity to possible 
non-SM physics. 

To incorporate the existing NNLO calculations, we 
proposed a new scheme that defines q^ independent coef- 
ficients, C7 and Cg, which are real in the SM. The Cy^g do 
not contain certain SM contributions (involving Ci_6,8)j 
which make the coefficients usually referred to in the lit- 
erature as Cy^ complex and dependent. We view C7, 
Cg, and Cio = Cio as the unknowns sensitive to physics 
beyond the SM to be extracted from data. 

Since precise measurements of the q^ dependences re- 
quire very large data sets, we studied how one can extract 
all Wilson coefficients obtainable from B Xgi^i^ 
from a few simple integrals of the Ht,l,a components 
of the decay rate. We concentrated on the low region 
(1 GeV^ < q^ < BGeV^) and found (see Figs. | and |) 
that splitting the total rate into transverse and longitu- 
dinal parts is a powerful tool to gain more information. 

The same angular decomposition in exclusive B — > 
K*t^t~ decay, together with the well-measured B 
K*j rate, also provides a determination of the Wil- 
son coefficients, without reliance on form factor models 
and without requiring a measurement of the zero of the 
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forward-backward asymmetry. In the heavy quark hmit, 
as a starting point of a systematic expansion, one can 
parameterize the form factors in the low region by 
just a few numbers [see Eq. ([l9[)]. Measuring the obscrv- 
ables in Eq. (|20|), one can extract from the data both 
the hadronic unknowns and the Wilson coefficients. A 
more complete understanding of the B K*i^£~ decay 
may be expected in the near future and will help to firm 
up the error estimates in such an analysis based on the 
B -> K*e+e- data. 

It is not known if a truly inclusive study of B — > 
Xdi^(~ will ever be feasible experimentally, but it may 
be possible to study the exclusive decay B p£^£~. 
The methods discussed in this paper are clearly applica- 
ble to these decays as well. 

With significantly more data, one may prefer to split 
the rate into more than two bins in the 1 GeV^ < < 
6 GeV^ region. One can then search for non-SM physics 
by fitting for complex Cy^g^io values, or allowing oppo- 
site chirality operators absent in the SM (model inde- 
pendently there is no preference which way to extend the 
parameter space). Given the good overall consistency of 
the SM, overconstraining determinations of the Wilson 
coefficients, such as that in Fig. ||b, may give the best 
sensitivity to new physics (similarly to the CKM fit). 

We did not include shape function effects jio) in our 
analysis for the inclusive decay (nor are they included in 
any other paper performing fits to extract short distance 
physics from the low region) . This is left for future 
work. Based on Ref . , we anticipate that if the B — > 
Xsj photon spectrum is used to understand the effect of 
the mx cut in B ^ Xsi^£~ , then the analysis considered 
in this paper will receive only modest corrections, leaving 
the general picture of how best to extract the Wilson 
coefficients unchanged. 



Acknowledgments 

We are grateful to Peter Cooper, Enrico Lunghi and 
especially Jeff Berryhill for helpful conversations. We 
thank Christoph Bobeth for providing us with his Math- 
ematica code for the running of the Wilson coefficients 
as given in Refs. [Ill M. Z.L. thanks the MIT GTP 



for its hospitality while part of this work was completed. 
This work was supported in part by the Director, Office 
of Science, Offices of High Energy and Nuclear Physics 
of the U.S. Department of Energy under the Contracts 
DE-AC02-05CH11231 (Z.L. and F.T.), and the cooper- 
ative research agreement DOE-FC02-94ER40818 (I.S.). 
I.S. was also supported in part by the DOE OJI program 
and by the Sloan Foundation. 



APPENDIX A: ANALYTICAL RESULTS 

Here we collect the explicit results for all and power 
corrections used in the main text. The Wilson coefficients 
Ci(^) refer to the operator basis of Refs. jl8[ except 
that wc keep Ct-iq in the traditional normalization [ p| . 
Their numerical values are collected in Table | below. 
(The corresponding coefficients in the normalization of 
Ref. [|l^ are as/(47r)C7_io.) Since the formally leading 
term in Cg(mb) is numerically small, it is often considered 
as 0(1) in the recent literature. In our case this is not an 
issue because we treat Cg as an unknown 0(1) parameter 
to be extracted from experiment. 

Throughout this paper we work in the IS" scheme [^, 
with rric^b always referring to the IS' masses ?tiJ"|. They 

pole 



are related to the pole masses, mP°^°, by (Cp — 4/3) 



,13 



pole 



An 



, (Al) 



where S^^{^) — {tt l2)as{^)CF is formally counted as 
0(1). By definition m}^ is /i independent. The ^ de- 
pendence reenters when the per turbative expansion on 
the right-hand side of Eq. (Al) is truncated. Switch- 
ing from the pole to the 15 scheme, one reexpands the 
perturbative expressions to a given order in as, count- 



ing 5^^ ^ 0(1), and using the same scale fj, in Eq. (Al) 



as everywhere else. All perturbative expressions below 
have been converted to the IS scheme. The expressions 
in the pole scheme given in the literature are recovered 



by setting 5^^ (fi) — and — m^°^° everywhere. 

In the following, repeated indices are summed from 1 
to 6. The coefficients Cy-io in Eq. ([Ill ) are defined as 



m 



Ci = C;7(At) To- + Ci[fi)Ki7 : In 



An 



Ki9 - TiQ In 



IS 



+ 0(«?): 



IS 



An 



(0) 



2 'J- 



+ 0(a2), 



(A2) 



where the higher order corrections are determined by the requirements that they vanish at /i = m^"^ and that the 
Ci are /i independent to a given order. For convenience we also defined Cg. It enters the Fi{q'^) (see below) at 0{as), 
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so we need it only at lowest order. The relevant entries of the anomalous dimension matrices are ||T^, [T^, [Tq) 



(0) 
lij 



/-4 


8 
3 





2 
9 





\ 


^(0) 


12 








4 
3 



















52 
3 





2 










40 


100 


4 


5 


' 'i9 


9 


9 


9 


6 













256 
3 





20 


/i9 







256 
9 


56 
9 


40 
9 


-V 



232 464 64 200 6464 11408\ 
243 '"sT' 8T'~243 ' SF ' 243~ j 
32 8 16 32 112 512\ 



27 ' 9 ' 9 ' 27 ' 9 ' 27 y 
2272 1952 6752 2192 84032 37856\ 



(A3) 



729 243 243 729 



243 



729 



The constants in Eqs. ( [A^ ) are included by convention, as in the usual definitions of , and are given by |l^ 



1 



20 80 



3 9 



«»7 = 0,0,--,--,-— ,-— , 0,0,1, --,20,-— , 



1 



10 



6 



Ki9 



0,0,0,1,^,1,, 



The functions Fi(q^) in Eq. ( pi] ) contain perturbative corrections arising from Oi-g.s- To 0{as), 



Aw 



[tjh{mc, s) + Ujh{l, s) + Wjh{Q, s)] 
J^Wi^m^J] , ^(s/4) 2, 



Q.(m), M (0) 
4tt 



47r 



G(/i)/f (m„s)+C8/f (s)' 



47r 9 



5^^(Ai)Q(M) 



l-s/(4?7i2) 'l-s/4 3 



-{t, + Ui + Wj) 



where s = /{^l^)'^ , ttT'c = m].^ /ra\^ , and 



1 , 6 , , 60 , 



7 2 32 
0,0, -38,-- 



0,0, 



32 

y 



(A5) 



(A6) 



The second line in Fg{q'^) arises from reexpanding the leading order contribution in the 15 scheme. The one- loop 
function h{rhc., s) encoding the four-quark contributions is ^] 



^ 8, _ , 4 2s 4,^ 
Mm„s) = --ln(m.)-- + ^--(l 

4 8 4 

/»(0,s)--gln(s) + — + -^7^, 



2ihi 



Ami 



and we have defined the function 



arctan 



< X < 1 , 



1-x' 

In (^+ Va!~~T) - i7r/2, x>l 



(A7) 



(A8) 



It satisfies J^{0) = 2/3 and - e) = sgn(e)(7r/2)yi -I- 0(e). 

The one-loop functions /g^'^"* (s) in Eq. ([A5|) containing the Og contributions are 
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4s 27r^ 

ln(s) + (22 + 60s 1583^ + 329s^ + ...)- (2 + 9s + 2As^ 50s^ + ..) 4i7r 

1 — s 3 



ln(s) + —(390 + 1480s + 3553s2 
1 — s 15 



27r^ 

- (4+ 15s + 36s2 + ...) — 



(A9) 



Their exact analytic expressions in terms of integrals can 
be found in Refs. ||l^, Here we expanded the terms 
not proportional to ln(s) for small s. Since these are small 



compared with the ln(s) term, the expressions above are 
accurate to better than 10~^ for s < 1 and 10~^ for 
s < 0.4. 
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The two-loop functions /^L's^ in Eq. ( [A5| ) contain the 
virtual 0{as) contributions from Oi_6 s-nd are known 
for i = 1 and 2 only. In Ref. they are given as 
expansions in s and rhc- Since the expressions for general 
rhc are quite lengthy, Ref. jl^ also quotes them at five 
fixed values for rhc, which we use together with a linear 
interpolation for intermediate values of rhc- 

Usually, the contributions from O1.2.8 that enter 

Cf^Q^{q^) are included via functions i^i^2 8^ [|4|; e.g., for 
the Oi contribution to Cn"'^' 



a 



incl 



47r 



(AlO) 



The i^i^2^g consist of terms containing powers of ln(/i/r7if,) 

and the functions f['^'fg- For example, switching to the 
pole scheme, the terms proportional to Q!s/(47r)Ci are 



(1) _ L,(o).,(o) 1„ JL 
rub 



+ j[f [tjh{7hc, s) + Ujh{l, s) + wjh{0, s)] I + /I 

\h{mc, s) 



(9) 



In 



nib 



256 ^^ fi 
243 ^m, 



2272 



729 



27 



hil,s) 



27 



h{0,s) 



(9) 



(All) 



for small s, we recover the result given in Ref. |Q. The 
same is true for all Note that in this way we au- 

tomatically include the ln(^/mf,) dependent terms of the 

analogous functions -Pjl'g ^ , which to our knowledge have 
not been calculated explicitly so far. 

The functions G7^Q{q'^) contain 0{l/m'^) corrections 
in Eq. ([Tl|), calculated for dF/dg^ and dAps/dq'^ in 
Ref. [12j. We found that their contribution can be in- 
cluded similarly to other four-quark operator contribu- 
tions in C™Q\q-^) [see Eq. (pi])] via the functions 



G9(g') 
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1 - 2s 



G7{q') 



6 / 1 — (7-^/(4™^) 



(A12) 



(In the operator basis used in Ref. Il2| the C2 — Ci/6 term 
in Eq. ( Al2| ) should be replaced by their C2, and we set 
this coefficient to unity in our numerical analysis.) These 
corrections diverge as {4m^—q^)~^/^ as q^ approaches the 
4m^ threshold. Following Ref. |l^, we imagine that this 
calculation makes sense for g2 < Sml « 6 GeV^ Even 
for q^ ^ 4to^, there is an infinite series of corrections 
suppressed only by increasing powers of Arub /to^ ■ 
The Luj{s) containing the 0{as) corrections in Eq. (|l 
can be extracted from Ref. p^. Defining 



In our scheme these ln(/i/m5) dependent terms are split 
up between _Fg(g^) and Ci, and are included in the defi- 
nitions (A2) and (|A5|). Expanding the term in brackets 



L{s) = 4Li2(s) + 21n(.s)ln(l-s) + 



2^2 



(A13) 



we find 



c^99(s) 



1 „n, , , l + 2s, ^ 5 + s(7-2s),^^ 19 + s l + 3sS^^{n) 

Lj^^{s)+L{s) + ln(l-s)-f — — i — r^ln(s' ' 



= 2 



2^ 
l + s(3-s) 



(1-s) 



s ■ ' (1 - s)2 ' 2(1 - s) 1-s 2 

T- / N l-s(9-2s), ,^ , 4s(5 - 2s) ^ . , 2(5 - 2s) , , „ 
2 Li2(s) + ln(l - s) - (l_^)2 Li2(V^) + 7^ ln(l - V^) 



■ S TT 



s(l-s) 



1 - s 



(1 



i cj^'*'^ (s) + L(s) + 3 ln( 1-s) 



1-s 2 ' 
s(l + 2s) 



ln(s) 



5 + s(12-s)^., 9-(2-Vi)%.,, ^ 
-4 — — ^ Li2(l + Li2(l - s) 



(1-s)^ 



5-|-s(47-4s) 3 + s(5i^(/i) 
2(1 - s) 1-s 2 

9 - (2 + Vi)2 

2 ' 



^ c.^^(s) + Lis) + 3 ln(l - s) - ln(s) 



(1 + %/^)^ 
16 + 5s(5-s) 3 -I- s ^^■^(/i) 



i= c.^^(s) + L(s) + ^ ln(l - s) + 2^^^^^ ln(s) 



6(1 -s) 1-s 2 ' 

21 - s 1 + 3s S^^ifi) 



12-(3-s)2 4_(i_^)2 
4 — — TTT^ Li2(l - Vs) 7-^ 7^7^ Li2(l - s) - 



(l-s)2 



(1 - V^sy 



2(1 - s) 1-s 2 
4 - (1 + V^)2 ^2 
(1 + V^)2 2 ' 



(A14) 



(A15) 
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- uj^^(s) + Lis) H ln(l - s) ) ln(s) - 

2 w -r V y T 2s ^ ^ 2(1 - s)2 ^ ^ 2(1 - s) 



1 + s 



3 + s(9-2s)^. , , l-s(22-s), . , 
^ Li2 (s) + — — ^ —!- ln(l - s) 



5(1 + s) 7r2 



2s(l - s) 
s + (3 - 2Vi)2 
2(l + yi)2 



2i±iiH_df)Li,(y7) 



1 



2(1 -s)2 3 
-u'^^s) + L{s) + i-t^ ln(l - s) 



1 - s 
s{7-s 



2 



(l-.s)2 



1 + lis 



(l-s)2 



Li2(l- VJ) 



1 + yi 

(1 - 



(l-,s 
■Li2(l-s) 



13- 3s 
1 - s 



ln(l - V^) 



2(1 - s) 1 - s 
(l + ^/i)2 2 • 



5^^(m) 



(A16) 



The functions (-u^{s) entering aj^^(s) contain all the de- 
pendence on -^/s, which cancels in the spectrum. All 
ln(/i/mt,) terms that usually appear in the functions 
ti;J^'^^(s) have been moved into Cy (along with the ap- 
propriate constant term contained in rni,/m}f). 

The xl{s) containing the 0{l/m1) corrections in 
Eq. (P) can be extracted from Ref. pi: 



x¥i'^ 
X'^{s 
X?i-s 



Ai +3A2 5 + 3s 



6 



1 - s 



2As 



s(4 - 3s) 



Ai + 3A2 3 + s(2 -I- 3s) 3 + s(4 - 3s) 

- 2A2 ■ 



6 (l-s)2 
Ai +3A2 3 + 13s 

6 1 - s 
Ai +3A2 3 + 5s 



(l-s)2 



2A; 



-2A2 



6 1 - s 
Ai +3A2 13 + 3s 



2A 



■(1-S)2' 

3-2s2 
s(4-3s) 



1 - s 
5-3s2 



6 

Ai + 3A2 

2 ' (1 - s)2 ' 

Ai + 3A2 3 + s(2 + 3s) 



(1-S)2 ' 

^ 5-)-3s(2-s) 



6 

Ai + 3A2 



(l-s)2 



(l-s)2 



-2A2 



(l-s)2 



(A17) 



APPENDIX B: NUMERICAL INPUTS 

In this Appendix we collect all of our numerical inputs. 
All values are taken from Ref. except where stated 
otherwise. To evaluate the Wilson coefficients we use 

mw = 80.403 GeV, 
sin^ Ow = 0.23122 , 

^poic ^ (171.4 ±2.1) GeV, 
Q^s(wz) — 0.1176 , 
= 80 GeV , 
/z^ = 120GeV. (Bl) 





fi = 2.35 GeV 


/i = 4.7 GeV 


/i = 9.4 GeV 


Qs(/i) 


0.2659 


0.2140 


0.1793 


ClM 


-0.4642 


-0.2880 


-0.1506 


C2M 


1.019 


1.007 


1.001 




-0.0096 


-0.0043 


-0.0017 


din) 


-0.1247 


-0.0795 


-0.0508 




0.00069 


0.00029 


0.00009 


C6(m) 


0.00205 


0.00081 


0.00026 


C8(m) 


-0.2012 


-0.1778 


-0.1598 


fnb{fj,) 


4.703 


4.120 


3.707 




-0.3637 


-0.3293 


-0.2982 


C7 


-0.2435 


-0.2611 


-0.2687 


C9(m) 


4.504 


4.209 


3.790 


Cg 


4.258 


4.207 


4.188 


Clo 


-4.175 


-4.175 


-4.175 



TABLE I: Values of the Wilson coefficients to 0{as) at dif- 
ferent low scales fi. 



Here, /ip * are the matching scales in the charm and top 
sector, respectively, and we use the same values as in 
Ref. iQ . For the top-quark mass we use the newest CDF 
and DO average |Q . The resulting values for the Wilson 
coefficients at 0{as) run down to the low scale and the 
corresponding values for the Ct according to Eq. ( [A^ ) are 
listed in Table ||. Note that the residual scale uncertainties 
of Cj and especially Cg are much smaller than those of 
C7,g(/i). We use a Mathematica code by Bobeth with 
the initial conditions and renormalization group running 
as given in Refs. |l9, EOl. For Cci{^) this requires the 



three-loop mixings calculated in Refs. |44 
In the decay rates we use 

aom("ih) 1/133 , 
|VtfcV,:.| =41.09x 10-3, 
ms = 5.279 GeV, 
tb = 1.584 ps, 
rriK' ^ 0.892 GeV, 



(4.70 ±0.04) GeV, 
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m^^ = (1.41 ±0.05) GeV, 
Ai = -0.27 GeV^ 
A, 0.12 GeV^ . 



(B2) 



We use the value of the electromagnetic coupling at the 
scale /z ~ rrif,, because for the total rate in this case the 
higher order electroweak corrections (which we neglect 
in our analysis) turn out to be numerically small, below 



The value of \VtbVt,\ 



the two percent level |20|, t 

taken from Ref. For ml^we take the naive average of 
Refs. which coincides with the PDG average [38|, 

and use the average of the errors quoted in Refs. |g^7 i5|. 
For m].^ we use the result of Ref. as quoted in the 
IS" scheme in Ref. pi. Finally, for Ai we take the value 



from Ref. 12711, and Ao 



(to 



;)/4 w 0.12GeV^ 
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